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Abstract 



> 

Oj^ I Existing calculations of heavy quark hadroproduction in perturba- 

^ ' tive QCD are either based on the approximate conventional zero-mass 

perturbative QCD theory or on next-to-leading order (NLO) fixed- 
flavor-number (FFN) scheme which is inadequate at high energies. 
We formulate this problem in the general mass variable-fiavor-number 
Pn\ scheme which incorporates initial/final state heavy quark parton distri- 

p • bution/fragmentation functions as well as exact mass dependence in the 

^ ■ hard cross-section. This formalism has the built-in feature of reducing 

to the FFN scheme near threshold, and to the conventional zero-mass 
parton picture in the very high energy limit. Making use of existing 
^ I calculations in NLO FFN scheme, we obtain more complete results on 

' bottom production in the general scheme to order both for current 

accelerator energies and for LHC. The scale dependence of the cross- 
section is reduced, and the magnitude is increased with respect to the 
NLO FFN results. It is shown that the bulk of the large NLO FFN 
contribution to the single heavy-quark inclusive cross-section is already 
contained in the (resummed) order "heavy flavor excitation" term in 
the general scheme. 
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I. INTRODUCTION 



The production of heavy quarks in high energy processes has become an increasingly 
important subject of study both theoretically and experimentally The theory of heavy 
quark production in perturbative Quantum Chromodynamics (PQCD) is more challenging 
than that of light parton (jet) production because of the new physics issues brought about by 
the additional heavy quark mass scale. The correct theory must properly take into account 
the changing role of the heavy quark over the full kinematic range of the relevant process 
from the threshold region (where the quark behaves like a typical "heavy particle") to the 
asymptotic region (where the same quark behaves effectively like a parton, similar to the 
well known light quarks {u,d,s}). Stimulated by significant recent experimental results on 
heavy quark production from HERA and the Tevatron, a number of theoretical methods 
have been advanced to improve existing QCD calculations of heavy quark production 
incorporating a dynamic role for the heavy quark parton. The purpose of this paper is 
to explain how the method of ACOT is applied to heavy-quark production in hadron- 
hadron collisions, to compare results of this approach to existing "NLO" calculations, and 
to demonstrate that it satisfies some important consistency conditions. 

Let us consider the production of a generic heavy quark, denoted by H, with non-zero 
mass rriH, in hadron-hadron collisions. We will define a quark as "heavy" if its mass is 
sufficiently larger than Aqcd that perturbative QCD is applicable at the scale rriH, i.e., that 
asimn) is small. Thus the c, b, and t quarks are regarded as heavy, as usual. Let Q be 
a typical large kinematic variable in the hard scattering process, in this case the px of the 
heavy quark (or the associated heavy flavor hadron). The details of the physics of the heavy 
quark production process will then depend sensitively on the relative size of the two scales 
tuh and Q. For simplicity, we will assume there is only one heavy quark H that we need 
to treat. Extending our treatment to the real world case of {c, 6, t} quarks with successive 
higher masses is straightforward. 

Conventional PQCD calculations involving heavy quarks consist of two contrasting ap- 
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proaches: the usual QCD parton formalism uses the zero-mass approximation {mu = 0) 
once the hard scale of the problem (say, Q) is greater than m^^, and treats H just like the 
other light partons on the other hand, most recent "NLO" heavy quark production 

calculations consider mn as a large parameter irrespective of the energy scale of the physical 
process, and treat H always as a heavy particle, never as a parton P|-pT|. We shall refer 
to the former as the zero-mass variable-flavor-number (ZM-VFN) scheme - since the active 
flavor number varies, depending on the energy scale fi ^ Q; and the latter as the fixed- 
flavor-number (FFN) scheme - since the parton flavor number is kept fixed, independent of 
Q. 

Each of these approaches can only be accurate in a limited energy range appropriate for 
the approximation involved: Q ^ rriu ioi the ZM-VFN scheme; and Q ~ ttlh for the FFN 
scheme. Nonetheless, both approaches have been used widely beyond their respective regions 
of natural applicability: on the one hand, NLO FFN calculations of c and b production are 
invoked from fixed-target to the highest collider energies ; and on the other hand, ZM-VFN 
calculations dominate practically all other standard model and new physics calculations, 
including most global PQCD analyses (which give rise to commonly used parton distributions 
0-^) and all popular Monte-Carlo event generators. The breakdown of the approximations 
beyond the original regions of applicability of these approaches can lead to unreliable results, 
and introduce large theoretical uncertainties. One possible sign of the latter is excessive 
dependence of theoretical predictions on the unphysical renormalization and factorization 
scale /i - which is well-known to be present in the NLO FFN calculation of hadro-production 
cross-section of c and b [lT2|,p^. 



With steadily improving experimental data on a variety of processes sensitive to the 
contribution of heavy quarks (including the direct measurement of heavy flavor produc- 
tion), it is imperative that the two diametrically opposite treatments of heavy quarks be 
reconciled in an unified framework which can also provide reliable theoretical predictions 
in the intermediate energy region, which in reality may well comprise most of the current 
experimentally relevant range for charm and bottom physics. This can be achieved in a 



general-mass variable-flavor-number (GM-VFN) scheme which retains the mn dependence 
at all energy scales, and which naturally reduces to the two conventional approaches in their 
respective region of validity [p!3|-p!5|J^ . The method is a development of the one devised by 
Collins, Wilczek, and Zee ||T^. The key point is that one can resum (and factor out) the 
mass singularities associated with the heavy quark mass ttih into H parton distribution and 
fragmentation functions without simultaneously taking the mn — > limit in the remaining 
infra-safe hard cross-section in the overall physical cross-section formula (as is routinely done 
in the ZM-VFN scheme[]). The resulting general formalism represents the natural extension 
of the familiar PQCD framework to include heavy quark partons both in the initial and 
final states of the hard scatterings which contribute to high energy processes in the SM and 
beyond. 

The principles and the practical application of this method were described in some de- 
tail for (heavy quark parton contribution to) Higgs production in Ref. and for lepto- 
production of heavy quarks in Ref. It has been applied to the analysis of charm 

production in neutrino scattering [IT^, and, recently, to a new global QCD analysis of par- 
ton distributions [|T^]. In the present paper, we apply this general formalism to heavy quark 
production in hadron-hadron colliders. For a concise summary of this formalism, and related 
recent developments, see Ref. for a systematic proof of the factorization theorem which 
provides the theoretical foundation of this formalism, see Ref. In recent literature on 
leptoproduction of heavy quarks, there have been two other formulations of the VFN scheme 
with non-zero heavy quark mass: an order scheme by Ref. Q and an order scheme by 
Ref. . For definiteness, we shall refer to our implementation of the general principles as the 
ACOT scheme. Although not unique, it represents in many ways the simplest and the most 
natural one in relation to the familiar MS ZM-VFN scheme See Ref. [501 for comments 



^The resummation of mass singularities into parton distributions and taking the zero-mass 
limit on the hard cross-section are usually performed simultaneously not as a matter of 
principle, but for the incidental reason that the mass singularities are most conveniently 
identified by using dimensional regularization in the zero-mass theory. 
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on some aspects of scheme dependence (in particular, mass-dependent or mass-independent 
evolution of the parton distributions), as well as on comparison to the scheme proposed by 
Ref. i. 

The following section presents details of our calculation scheme as applied to hadropro- 
duction and describes the physical origin of the various terms which appear in the formalism. 



Sec. m contains specific information on how the various ingredients of the formalism are 
calculated. This is followed by numerical comparisons of the new calculation with existing 
FFN scheme results on the inclusive differential pt cross section for bottom production at 
current accelerator and LHC energies. In the concluding section, we discuss what remains 
to be done for a full understanding of heavy quark production in PQCD. 

II. HADRO-PRODUCTION OF H IN THE GENERAL MASS FORMALISM 

Consider the hadro-production of a generic heavy quark H: 

A + B — >H + X (1) 
where A and B are hadrons, and X includes H along with all other summed-over final state 



particles. In the ACOT formalism for heavy quark production developed in Refs. |13i|2|,p!9 
the inclusive cross section for observing H with a given momentum p at high energies is 
given by a factorization formula of the same form as in the familiar zero-mass QCD parton 
formalism: 

^ff = E ffs ® ® rff , (2) 

a,b,c 

where {A, B} denote the initial state hadrons, {a, b} the initial state partons, /^^ the asso- 
ciated parton distribution functions in the combination /^^ = /^/^, d'^b ^^e perturbatively 
calculable infra-red safe hard-scattering cross section for {a, b} {c, r} which is free of large 
logarithms of mn over the full energy range, and the fragmentation functions for finding 
H in c. All active partons are included in the summation over {a, b} and {c}, including H 
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provided the renormalization and factorization scale /i is larger than mn-^ 

Since we shall compare our results to those of the FFN scheme, it is necessary to draw the 
distinction between the heavy quark H and the associated light partons, we shall denote the 
latter collectively by I for simplicity of notation. By definition, I = {g,q} where g is the gluon 
and q denotes the light quarks in the sense of the FFN scheme: q = {u, d, s} / {u, d, s, c} for 
charm/bottom production respectively. The number of light quark partons will be denoted 
by ni. 

As mentioned in the introduction, an important distinguishing feature of this general 
formalism from the familiar ZM-VFN parton approach is that, after subtraction of mass- 
singularities, the full niH dependence in the perturbatively calculated hard scattering co- 
efficients (T^^^ is retained. This allows the theory to maintain accuracy and reproduce the 
FFN scheme results in the threshold region, as required by physical considerations. On 
the other hand, our parton densities (and fragmentation functions) are defined in the MS 
scheme, hence they satisfy mass-independent evolution equations — the same as in the 
usual zero-mass formalism. Considerable simplification then results in the implementation 
of this scheme since the well-established NLO evolution kernels and evolution programs can 
be directly used. It should be noted, however, that the parton densities do have imphcit 
dependence on quark masses: to the leading power in A/mn, this dependence is generated 
by the matching conditions at yU = mn between the parton densities below and above the 
threshold]^ — they are defined in each of the two regions by the respective renormalization 
scheme adopted for that region [ p!3| , p!^ . 

The first few terms in the perturbative expansion of the production cross section, Eq. |^, 



^For simplicity, we shall use the symbol fi to represent collectively the renormalization scale 
as well as the factorization scales for parton distributions and for fragmentation functions. 



^ Note that equivalent matching conditions can be derived for any /i which is of order tuh- 
Also, possible non-perturbative heavy quarks as opposed to "radiatively generated" 
ones (assumed in subsequent discussions), can be incorporated in the general scheme by 
allowing for a nonzero heavy quark density in the below-threshold part of the scheme. 
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are schematically 



fVr, (5?) 
JAB ^ 












+ffB^ 
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+ flB 










: HCl 


+ flB^ 










: GF1,2 












: HH2 



where repeated indices I, I', h, and I2 are summed over all light flavors {q and g) and the 
pre-superscript n on "o" denotes the formal order in the expansion of a^. 

The physical interpretation of the various terms, labeled by the abbreviations in the last 
column of this equation, can be made apparent by the generic diagrams shown in Fig. |l|. 
They are: 

(i) HCO: the order heavy flavor creation process I V ^ H H followed by heavy-quark 
fragmentation; 

(ii) HEl: the order a] heavy flavor excitation process I H ^ I H followed by heavy-quark 
fragmentation; 

(iii) HCl: order virtual and real corrections to HCO, / /' HHg, with fragmentation; 

(iv) GF1,2: order (light- and heavy-) parton-gluon scattering I g —>■ I g and H g ^ H g, 
followed by gluon fragmentatior^ into the heavy quark if; and 

(v) HH2: the order a1 HH scattering HH{H) HH{H), with fragmentation. 



4por generality, we have written the summation (over /) to include all light-parton frag- 
mentation into H. In fact, we expect gluon fragmentation to dominate. Light-quark frag- 
mentation is suppressed because of the quark structure of the evolution equation. In the 
flxed order calculations which we will do to compare schemes, fragmentation of light quarks 
q and heavy antiquarks H into H are both at least of order a^. 
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HCO 



HE1 



HC1 





GF1 



GF2 



HH2 



legends 

- gluon / light-quarks 

— heavy quark 



I I hard cross-section 
parton distr. fn. 
fragmentation fn. 



FIG. 1. Graphical representation of the leading terms in the factorization formula which 
correspond to the various production mechanisms. The initial state hadron line for the parton 
distributions are uniformly suppressed. 

For the diagrams in Fig. |l] we use a dashed line for all light partons / (representing 
collectively gluons and light quarks) to distinguish them from the heavy quark if, represented 
by a solid line. The square blocks in these diagrams represent the hard cross sections "■a, 
obtained from all Feynman diagrams of the same order with the external lines indicated. 
The round blobs represent parton distribution and fragmentation functions. The initial state 
hadron lines for the parton distribution factors are suppressed in all the diagrams. For the 
differential cross section da/dp^, the four terms involving 2-^2 hard scattering processes 
are tree-level processes. Beyond tree level, we have included only the most important order 
al term — HCl — which contains 2^3 real and 2 — > 2 virtual corrections to the leading 
HCO cross-section. This term corresponds to the next-to-leading (NLO) contribution in the 
FFN scheme; it is known to be large and it is the only one that has been calculated to this 
order so far (see below). We will discuss the relative sizes of the various terms in the next 
paragraph, and comment on the possible significance of other terms not included here in 
later parts of the paper. 

In conventional applications of PQCD with light partons, it is common to distinguish 
the order of terms in a cross section by treating the (non-perturbative) parton distributions 
in as all being of order one, and then counting powers of as in the hard scattering. 
With the inclusion of heavy quark partons, one may expect heavy quark distribution and 
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fragmentation functions to be of order in the energy region not too far above threshold, 
assuming they are purely radiatively generated.^ We know that at large x the valence u 
and d quarks dominate all the other partons; while at small x, it is the gluon density that 
dominates. For instance, in the latter region, we compare in Fig. ^ the relative sizes of the 
various sea partons at two energy scales to the gluon (scaled down by a factor of 5). In all 
cases it appears safe to regard the fully evolved heavy quark density fj^ to be of effective 
order with respect to the dominant parton density (gluon or valence quark), i.e.. 



fA ~ 0{as 



(4) 



Note, however, that the evolved charm density comes within a factor of 2 of the other sea 
quark densities at small x. 




H = 10 GeV 




H = 100 GeV 



0.001 0.003 



0.01 Y 0.03 



0.001 0.003 



0.01 Y 0.03 



FIG. 2. PDFs vs. X for ^ = 10 and 100 GeV. The gluon PDF is scaled by 1/5. 



The same argument leads to the following expectations for the fragmentation functions 
(note, we confine ourselves to the production of the heavy flavor parton H, rather than the 
associated hadron): 



d" 



5{l-z) 

0(«s) 



0{as 



(5) 



^As mentioned in footnote |^, our formalism can accommodate non-perturbative, i.e. non- 
radiative or "intrinsic" heavy quarks. However, we shall not get into that possibility in this 
paper. 
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These order-of-magnitude estimates are, of course, expected to become inapplicable at very 
large scales. The gluon-to-heavy-quark fragmentation function, in particular, can become 
substantial at large scales. 

Nevertheless, let us use these estimates as the first guide to orders of magnitudes of the 
terms in Eq. |^. The first term (HCO) is of order a^; the HEl and GFl terms, and the HCl 
term are of effective order a^; and the GF2 and the HH2 terms are of effective order a^. 
(This naive counting of effective order explains the choice of numerical suffixes in the labels 
for these terms.) The actual relative numerical importance of the various terms will also 
depend on other considerations such as large color factors or dynamical effects (e.g. spin-1 
t-channel exchange contribution at high energies). This is well-known for the perturbative 
expansion of various processes, including heavy quark production calculated in the FFN 
scheme, where the order NLO term (corresponding to HCl without heavy quark mass 
subtractions) is actually larger than the order LO term (HCO). We will examine in 
detail the numerical significance of the various terms in Sec. 0. The results are revealing, 
since they show simple features which are not accessible from the conventional FFN scheme 
perspective. 

To specify fully the calculation in our scheme, we need to specify the perturbative hard 
cross sections a'^^ . Following Ref. |^, we obtain these by applying the factorization formula, 
Eq. H, to cross-sections involving partonic beams, that is to a'^^ (without the caret) calculated 
from the relevant Feynman diagrams. Then we solve for the hard cross sections a'^'^ order 
by order in a^. Specifically, in a given order a", the factorization theorem implies that the 
partonic cross section has the form 

vf/ = E "^^J® (6) 

where the sum is over the values of rii, n2, and 77,3 that give the correct order, i.e., n = 
ni+n2 + n^. On the right-hand side, this equation differs from Eq. ^ in two respects. First, 
the parton distributions are relative to an on-shell parton target, instead of a hadron target. 
Secondly, we have expanded the parton-level distribution and fragmentation functions in 
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powers of ag, with, for example, "^/^ denoting the term of a"^ in the distribution function. 
These coefficients are calculated order by order in a^, but are not infra-red safe; they are 
used only in the intermediate steps toward the derivation of a^'^. 

The most general method to obtain the hard scattering coefficients is directly from Eq. 
^. All of the full partonic cross sections, the left-hand side, and the parton densities and 
fragmentation functions at the partonic level, on the right-hand side, can be computed from 
definite Feynman rules. From them, Eq. ^ gives the hard scattering coefficients. However, 
at the order level to which we are working, it is convenient to carry out this calculation 
in two steps; the relation to the coefficient functions of the other schemes in the literature 
is then readily obtained. 

First, the unsubtracted cross-sections are calculated from all relevant Feynman diagrams 
with the use of dimensional regularization (with the familiar parameters e and /i), with 
all light quark masses set to zero, and with the heavy quark mass non-zero. After ultra- 
violet renormalization (with MS counter terms) and cancellation of infra-red divergences 
(after combining real and virtual diagrams), the cross-section formula will depend on the 
renormalized mass rriH and the unphysical parameters (e and /i), in addition to the kinematic 
variables. The collinear singularities associated with the 1/e poles can be factorized into 
a kind of distribution function of light partons in a light parton f^{e,fi) using the MS 
convention. Since the cross-section is inclusive with respect to the light partons, we need no 
light parton fragmentation functions. Thus,0 

V2f(mH,e,/x)= J2 "V:J(e,/x)® "^af,^(mH,/x), (7) 

ni+n2='n 

where /^^ represents the product of two singular light parton distributions (cf. Eq. ||), 
and all non-essential (kinematical and convolution) variables have been suppressed. This 
concludes the first step of the construction, where the set of intermediate finite cross-sections 



^For simplicity, we have used the symbol fi to collectively represent both the factorization 
scale and the renormalization scale. 
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"'a^fj^ {niH, ^) are obtained by systematically subtracting the coUinear singularities. 

At first sight, this appears to imply that the sub-set of "cr^^(m//, /i) with {a,b} = 
{I, I'} (light partons) are the same as the conventional FFN scheme cross-sections. This 
is in fact true at the lowest non-trivial order, which is all that we will consider in this 
paper. However, at sufficiently high order extra heavy-quark loops come in, and these 
are renormalized differently in the FFN scheme and in the MS scheme that we use when 
/i > rriH- Hence the singularities to be subtracted differ by some kind of renormalization- 
group transformation. In this paper, we do not treat these higher order graphs. 

The second stage of our calculation starts from the observation that although they are 
finite, the cross-sections "a^"^(m//, /i) contain logarithmic mass-singularities, i.e. powers of 
In (m/^//i), in the mu/ ^ limit. The second step of the derivation consists of factoring 
these singularities out to arrive at the fully infra-red safe hard cross-sections "a^J,*^ of Eq. |^. 
Explicitly, 

Here the logarithmically singular terms in the rriH — > limit are factored into / and d. 
These are partonic level parton densities and fragmentation functions with subtractions 
made to remove the singularities associated with light partons. These subtractions exactly 
correspond to the subtractions used to obtain a from a. The remaining infra-red safe mn 
dependence is kept in a (in contrast to the conventional approach, where mn is set to zero). 

As with all calculations of hard scattering coefficients, there is a freedom to choose 
exactly how to define the parton densities and fragmentation functions. The choice defines 
the factorization scheme, cf. [^], and thus determines how much of the finite mn dependence 



is included in / and d. 



We use the ACOT scheme [0,0,|T^]. In this scheme, the parton densities and fragmenta- 
tion functions in the region fi > mn are determined by the requirement that all ultra-violet 
divergences are renormalized by the MS scheme — both the ultra-violet divergences in the 
renormalization of the interactions of QCD and the ultra-violet divergences that make finite 
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the parton densities (in hadrons). The hard scattering coefficients a in Eq. ^ are then well- 
defined. In the first stage of our calculation, the factorization of light parton singularities, 
we define the intermediate coefficients a in Eq. ^ by subtraction of coUinear singularities in 
the MS scheme. 

Subtracting the terms which contain / and d with non-trivial mass singularities from 
"'a^^{mH, fJ'), we obtain the fully infra-red safe hard cross-sections "'a'^'^rriH, fi) which we 
need in Eq. ^. The relevant non-vanishing perturbative parton distributions "/ up to order 
as, are: 



Vf(x,/i) = ^ In(^) P^^nix) (9) 
V|(x,/i) = ^ In (4) 6(1 -X) 

where a,b = {g,H}, Pa^b{x) is the usual first order splitting function. The only nonzero 
terms at order a<j in Eq. |^ come from contributions with a heavy quark loop. All one-loop 
corrections that involve light partons are zero, because of the well-known cancellation of IR 
and UV singularities. For the fragmentation functions, we have [22]: 



'dUz) = 5(1 - z) 

' ^ ' l\n{l-x) - 1^ (10) 



'dUz) = ^ 



1-x \ Mfj 



We can now substitute Eqs. || and |T0| into Eq. || to solve for ^o'^l- To order a^, only ^ f\ 
and ^c?!! are needed and we obtain the trivial identities, 



l-HH _ 2-HH 

"i I' — I' 

2~IH _ 2f.lH 



2~HH _ 2-HH 



^ HH ^ HH 



(11) 



which refiect the fact that the order aj. cross sections are given by tree graphs: i.e. the ex's 
are infra-red safe, in the mu — > limit. Hence they do not need any subtraction. Following 
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the same procedure to order a^, we obtain schematically:^] 

V 3~HHX _ 3±.HHX 

C^n'(real+virtual) ~ 



+ 
+ 



//' ® (^iH + Ji ® (^Hi' + ® d, (12) 



a 1 1, ® djj 



For convenience, we have used a single equation to cover the various possible initial states. 
Not all terms on the right-hand side are applicable to all cases: for gluon-gluon scattering, 
{/ /'} = {gg}i all terms are present; for {/ /'} = {qq} , only the 1st, 3rd and 5th terms 
contribute; and for {//'} = {gq{q)} , only the 1st, 4th and 5th terms contribute. 

This equation can easily be inverted to obtain the order hard cross section in terms 
of the calculated finite intermediate cross-sections (with non-zero mass mH) and various 
subtraction terms which remove the mass singularities associated with the heavy quark 
degree of freedom: 



^X (Til 



3-HHX 



V 3;^HHX 

^X 0//'{rcal+virtual) 




: HCl-FFN 


1 M ^ 2~HH 1 fg' . 




: Fin-Sub 






: HFl-Sub 






: HEl-Sub 






: GFl-Sub 



(13) 



Here, we have replaced all tree-level in Eq. |TT] by the corresponding ^5" because they 
are the same, cf. Eq. |ll|. The content of the terms on the right-hand-side of this equation, 
labeled by the abbreviations in the last column, can be seen more easily from the diagrams 
in Fig. ^ where, for clarity, the four possible initial/final state channels are separately shown. 
The terms in Eq. O and Fig. 13 are: 



^It should be observed that the only contribution to ^/|, in Eq. |^, corresponds to a gluon 
self-energy graph on the external gluon line. Therefore when obtaining a from a in Eq. |l^, 
the effect of the ^/^ terms is simply to cancel external line corrections. 
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& 2^3 real 



Fin HF1 HE1 GF1 

mass subtractions (2^2 collinear configuration) 












O intermediate partonic X-sec. q's 

ligint quark 

Ineavy quark 

FIG. 3. Graphical representation of the terms in the right-hand-side of Eq. 13 



'TTSTSTTB^ glUOH 



perturbative parton 
distri./ fragm. fn's 

lieavy quark loop 
self-energy corrections 



Collinear 

singularities due to light partons have already been subtracted. The vertices represented by a dot 
are the heavy quark parts of the perturbative distribution and fragmentation functions, as in Eqs. 
|g| and ITOI. 
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• HCl-FFN : the usual order FFN scheme result, due to contributions of the NLO- 
HC (virtual and real) diagrams, with infra-red and collinear singularities associated 
with light partons cancelled/subtracted in the conventional way (in the MS scheme); 

• Fln-SUB : the correction to the order gluon-gluon HC process due to the difference 
in the definition of the gluon distribution between the below-threshold {uf = ni) 
scheme and the above-threshold (n/ = ri; + 1) scheme; 

• HE-SUB : the subtraction of large logarithms of the heavy mass contained in ^/^, 
cf. Eq. due to a flavor-excitation configuration; 

• HF-SUB (GF-SUB) : the subtraction of large logarithms of the heavy mass in the final 
state heavy-quark (gluon) fragmentation residing in ^d^g, cf. Eq. |10[ 

After these subtractions, ^d'fi,^^ is free from all large logarithms associated with potential 
mass singularities, i.e. it is infra-red safe with respect to mnQ 0. When this result is used 
in Eq. |^, the hadronic cross section is well-behaved in the high energy limit, in contrast to 
the FFN scheme result which would diverge because of the large logarithms. In fact, in this 
limit Eq. § reduces to the usual zero-mass MS parton formula with the H quark counted 
just like the other light partons. This is, of course, the correct limit at high energies. 



Some insight can be gained by explicitly substituting Eqs. |TT|, |T3| in Eq. ^ and obtaining 
the hadronic cross section in terms of the intermediate partonic cross-sections {^a,^ a} (which 
are all finite for non-zero itlh) along with the necessary subtraction terms which represent the 
overlap between the two sets of cross-sections and which remove the potentially dangerous 
mass singularities. The full set of terms are most clearly displayed in diagrammatic form, as 
shown in Fig. ^. The first column lists contributions from all the 2 —>■ 2 tree-level diagrams 
summarized in Eq. ^ and Fig. ^ plus 1-loop virtual corrections from Eq. ^ the last column 
lists contributions from all the 2^3 terms contained in Eq. |T^; and the middle column 
contains all the relevant subtraction terms. These terms can most easily be obtained by 
substituting the terms shown in Fig. |^ in those in Fig. |I], summing / over {g,q}. With the 
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2^2 



heavy quark mass subtr. 



+ 



HC0,1 



HE1 



GF1 



GF2 



HH2 




FIG. 4. Graphical representation of the physical cross-section, Eq. |^ and Fig. ||, written 
in terms of the intermediate partonic cross-sections "cj's and the attendant heavy quark mass 
subtraction terms which represent the overlap between the 2^2 and 2^3 cross-sections. Initial 
state parton distribution function factors are uniformly suppressed. 
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exception of uniformly suppressing the initial state parton distribution factors (represented 
by dark blobs in Fig. |I|), these diagrams contain all the ingredients needed to write down the 
full formula for the cross-section. In between the last two columns, we have also drawn lines 
to indicate the origins of the various subtraction terms in the 2^3 diagrams, according to 
Fig. |. 

In this way of organizing the results, the order as subtraction terms in Fig. ^ are shown 
next to the relevant HE/GF/HF contributions in the same row, making explicit the phys- 
ical origin of the subtractions: these terms containing large-logarithm (present in the un- 
regulated FFN scheme calculations) also represent the low-order components of the QCD 
evolved parton distribution/fragmentation functions. As an example, consider the first two 
terms in row one of HEl: 

(2 2d) - (m^-subtraction)l^^;^''' = /1(/J - /I ® 7f ) ® '^f/ ® 4 (14) 

Both from this formula and from the corresponding graphs, one can see that the subtraction 
term containing ® ^ f^^ represents that part of the NLO-FFN contribution which is already 
included in the (fully evolved) parton distribution f^. The latter, of course, represents 
the result of resumming all powers of In (/i^/m|^), hence contains important physics not 
included in the NLO-FFN calculation, in addition to being well-behaved as ^'^ /m^j becomes 
large — according to the renormalization group equation. Similar comments apply to the 
other 2— *2 a terms and the corresponding subtraction. The GFl terms will be proportional 
to {d^ — ^d^ ®d^) where d^^ are the fully evolved and ^d^ the perturbative fragmentation 
functions. 

The following (terms in) parton distribution and fragmentation functions ^/^, ^/|, 
and d^ all vanish at the 'threshold', fi — ttih by calculation (cf. Eqs. |^JT^ and Refs. |T^J^). 
In addition, in the threshold region, 

d^ S{1 - z) + 0{as) 
fA-fA® 7f - 0{al) (15) 
(if - ^Jf ®4~0(«^) 



Hence, the differences between the first two columns of the terms represented in Fig. ^ — 
(2^2 5")— heavy quark mass subtraction — vanish even faster than the individual terms 
approaching the threshold. Using these results, we obtain, in this limit 

^HX threshold ^ 2~HH . LO-FFN 

+ fls ® ( ^^z^'Vtuai) + ^^/f'lLi)) : NLO-FFN (16) 
(full NLO - FFN scheme) 

That is, the hadronic cross section in this formalism reduces to the full flavor creation (HC) 
result of the FFN scheme to order a^. In this region, there is effectively only one large 
momentum scale {pr ~ ttih)] and the FFN scheme is well suited to represent the correct 
physics. 

We see, therefore, the ACOT formalism provides a natural generalization of the familiar 
light-parton pQCD to the case including quarks with non-zero mass which contains the 
right physics over the entire energy range. At high energies, Eq. |] gives the most natural 
description of the underlying physics. The mass subtraction terms appearing in Eq. |l^ for 
the hard cross section correspond to the poles arising from coUinear singularities 
in the mn = QCD parton formalism which are usually removed by MS regularization.0 
On the other hand, at energy scales comparable to the quark mass mn, consistency with 
the physically sensible NLO-FFN calculation in the threshold region is also guaranteed 
by keeping the full rriH dependence in the partonic cross sections appearing in Fig. ^, as 
described in the previous paragraph. 



^Consistency of this generalized formalism with the usual MS scheme at high energies is 
ensured by adopting the precise definitions of as(/i) and f^^, in the renormalization 
scheme of Ref. [|l3i|T6| . As mentioned earlier, in this scheme, the parton distributions and 
fragmentation functions satisfy the familiar MS evolution equations above the respective 
thresholds. 
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III. CALCULATIONS 



Our numerical calculations are carried out using Eq. 0, with hard cross sections given 



by Eqs. |TT] and |T3|. (The equations are presented graphically in Fig. ^.) We briefly describe 
how the various quantities on the right-hand-side of the equation are obtained. 

For the parton distributions fXi^yQ)^ we use the CTEQ3M set |^ for definiteness. The 
CTEQ parton distribution sets, in general, are evolved according to the ACOT scheme 
described in Ref. 0,|IB,R,r^ which is the one we use in the heavy quark production theory. 
This is necessary for obtaining consistent results (a point which has not been entirely obvious 
to all users of the scheme). The reason is: the expected compensation in the threshold 
region between the QCD-evolved Q) and the perturbatively generated subtraction 

fAi^jQ) SUB = Ia^^^Q) ® ^f^i^^Q) (cf- Eq- HID will not take place unless the choice 
of evolution scheme and the choice of the location of the heavy quark threshold match 
properly.0 This is illustrated in Fig. ^ which compares /i) and /^(x, /x) sub for charm 

and bottom as a function of /i at x = 0.1 and x = 0.01. Each curve individually vanishes at 
the threshold fi = mn required in this general scheme. As fi increases, both grow at a 
rapid rate because the evolution is driven by the large gluon distribution (through the mass- 
independent splitting kernel); however, the difference of the two, which determines the actual 
correction to the main contribution in this region (due to the HC process), grows slowly as 
one would expect on physical grounds. A failure to ensure the proper compensation between 
these terms due to a mismatch of schemes or the location of the threshold can lead to quite 
unphysical results because then the difference would be of the same order of magnitude as 



^As an example, MRS distributions use fi = 2mH as the threshold for evolving f^{x,fi), 
from zero in contrast to fi = rriH which is required by our renormalization scheme. Using 
MRS distributions in this framework represents a mismatch of schemes, and will lead to 
unphysical results as described below. If a matching point other than fi = mn is used, then 
the heavy quark distribution must start from a non-zero value. 
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FIG. 5. Comparison of the evolved PDFs, /^(x,/x) (labeled PDF), and perturbative PDFs, 
^f^[x,fi) (labeled SUB), as a function of the renormalization scale /i for charm at x = 0.1 (a) 
and X = 0.01 (b), and for bottom at x = 0.1 (c) and x = 0.01 (d). This shows the compensation 
between fully evolved heavy quark parton distribution and the first order perturbative contribution 
(which is the only part contained in the FFN scheme calculation). 



In calculating the partonic cross sections which appear in Eqs. |TT] and |T3|, the formula for 
the LO-FFN cross sections ^d'^g^ and ^d'g^ are well-known. For the NLO-FFN cross section, 
both virtual ^5"^^ and real ^ct;^^' ", we used the Fortran codes from Ref. |11] and Ref. [^. 

To calculate the gluon and heavy-quark fragmentation functions appearing in Eq. ^, we 
need also the QCD evolved fragmentation functions d^^^. As mentioned earlier, for the 
purpose of this paper, we restrict ourselves to the if-quark production cross section; hence 



^^f^{x, fi) and f^{x, /i) sub are not expected to cancel at large /i; the latter, in fact contains 
the divergent ln{fi/mH) factor in that limit. In that region, this divergent subtraction term 
plays the other important role of cancelling the corresponding large logarithm in the NLO- 
FC term to render the latter infra-red safe, as shown in Eq. |1^. 
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dgff are parton to parton fragmentation functions. To obtain heavy flavor hadron production 
cross sections, it is necessary to perform an additional convolution with the appropriate 
hadronic fragmentation functions. The QCD evolved fragmentation functions d^jj{z,fi) 
are generated by solving numerically the QCD evolution equation in our scheme, using as 
input at fiQ = irtH the perturbative formula, cf. Ref. The comments about proper 

compensation between parton distributions made above also apply to the fragmentation 
functions d^^ and ^d^^. These features have been examined in detail during our calculation. 

IV. RESULTS 

We now present typical results for h quark production cross section ^^^^ \y=Q vs. pt and 
vs. the QCD scale parameter fi at collider energies. For simplicity, as is customary in the 
literature, we use a single scale parameter fi to represent the renormalization scale, the fac- 
torization scale for the parton distributions, and the factorization scale for the fragmentation 
functions. In principle, these could be chosen as independent; the hard cross sections would 
then depend on all three scales. As a rule, we shall express the scale /i as multiples of the 
natural physical scale = ^Jpf + m^, although, again, other choices could also be con- 
sidered. Except for explicit discussions concerning the //-dependence of the cross sections, 
our default choice of scale is //^ = In order to be able to clearly discern the various 

contributions to the steeply falling function ^p^^y iPt), we shall in general use the scaled cross 
section d{pt) = Pt dp^dy \y=o when examining its pt behavior .0 We concentrate mostly on 
b-production at the Tevatron for definiteness. 



^^We also note, in evaluating the right-hand-sides of Eq. ^ (cf. also Fig. |), we have uniformly 
omitted the last convolution with d^ for simplicity of calculation. The numerical effect of 
including this factor is relatively small. 
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A. Inclusive pt distribution and comparison of heavy parton picture in the FFN 

scheme 



Fig. |]a shows o"(pt) vs. pt for h production at 1800 GeV, including the individual terms 
on the right-hand-side of Eq. Over the range < pt < 100 GeV, the two largest 
contributing terms are "leading order" (2 2, tree-level) HEl {o-gH^gu) and HCO {o^gg^Hs)] 
the other tree-level terms and the "next-to-leading" term, HCl {cFgg^gHH^ mass-subtracted), 
constitute about 10 — 25% of the cross section, depending on the value of pf The left-hand- 
side plot shows the contributions from the individual terms; and the right-hand-side plot 
compares the relative sizes of the combined LO-VFN (i.e. tree-level) contribution, the NLO- 
VFN contribution and the total cross-section TOT-VFN when the calculation is organized 
in this scheme. 




FIG. 6. Contributions to the scaled cross section Pt da / dpf /dy\y=Q (nbGeV^) vs. pt for b 
production at 1800 GeV with /u = Mt/V^ organized in the ACOT formalism, a) The curves 
correspond to the separate terms of Eq. ^. b) The curves are leading-order (a^) (LO-VFN), 
next-to-leading-order (af) (NLO-VFN), and the total result (TOT-VFN). 



Two interesting features are worth noting. First, the LO-VFN contributions (tree pro- 
cesses) give a reasonable approximation to the full cross section, the NLO-VFN correction 
is relatively small. (This is in sharp contrast to the situation in the familiar FFN scheme 



^^In terms of Fig. ^ the mass-subtraction terms are combined with the associated order 
HCl terms (from which they originate) to yield the infra-red safe hard cross-sections a. 
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where the NLO-FFN term is bigger than the LO-FFN one. Cf. Fig. ^ and discussions be- 
low.) This is, of course, an encouraging result, suggesting that the heavy quark parton 
picture represents an efficient way to organize the perturbative QCD series. Secondly, the 
HEl contribution is comparable to, and even somewhat larger than, the HC0 one - in spite 
of the smaller heavy quark parton distribution in the initial state compared to the gluon 
distribution. Closer examination reveals that two effects contribute to this non-apparent re- 
sult: a larger color factor for the HEl process, and the presence of t-channel gluon exchange 
diagrams which is absent in the HC0 process.0 




FIG. 7. Contributions to the scaled cross section Pt da / dp^ / dy\y=Q (nbGeV^) vs. pt for b 
production at 1800 GeV with = Mt/V^ organized according to the FFN scheme, a) The 
curves correspond to the order (LO-FFN) and order (NLO-FFN) heavy-flavor creation (HC) 
contributions (without heavy mass subtractions). Also shown are the corrections due to the HE, 
GFl, GF2, and HH2 terms with associated subtractions as given in Fig. ^. The last three are 
numerically negligible and appear at the bottom of the plot unlabelled. Note that, (i) the NLO 
term is two times larger then the LO one; (ii) the contributions from are small and unlabeled, 
b) LO-FFN and NLO-FFN contributions along with the total result (TOT-FFN). Cf. Fig. || for 
comparison with the ACOT scheme case. 



^^The precise values of the HEl contribution are somewhat sensitive to the choice of fac- 
torization scheme and scale, especially close to the threshold region, as will be shown below. 
However, one of the important features of our formalism is that any scheme and scale de- 
pendence in HEl will be closely matched by changes in the HCl contribution (through the 
corresponding subtraction term), so that the combined inclusive cross section remains rel- 
atively stable. Cf. discussions in the previous section about the matching of evolved and 
perturbative parton distributions. For the current discussion, we adopt /i = Mt/\/2 as a 
central choice, given commonly used ranges of fi such as [Mt/2,2Mt] and [Mr/4,Mr]. 
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It is useful to compare the above situation with the same results organized in a way more 
familiar from the conventional FFN scheme point of view. For this purpose, one begins with 
the intermediate cross-sections a for heavy flavor creation processes only, HCO {d^gg^Hil) ^i-nd 
HCl {o'gg^gHH, uo masS" sub t r act lou ) . Corrections to the FFN scheme calculations in the full 
scheme then consist of the remaining terms on the right-hand-side of the cross section formula 
depicted in Fig. ^, which are now most naturally organized with the mass-subtraction terms 
combined with the corresponding 2 — 2 cross-sections in the same row. Fig. ^,b show <T{pt) 
vs. pt in the same format as in the previous plot but with individual contributions organized 
in this way.[^ The largest term is now the "NLO" HCl {^gg^gHi}) followed by the "LO" 
HCO of the conventional FFN scheme. The fact that the NLO (order a^) term cthci is much 
larger than the LO term (order a^) anco ^ the "K-factor" is typically of the order ~ 2.5 - 
is disturbing from the perturbation theory point of view, as has been known since the order 
al calculations were first done. On the other hand, we see from Fig. ^a, the corrections 
to the FFN scheme terms, consisting of the other terms in Fig. ^, are positive but not very 
large - again of the order of ^ 20%. This means that the effects of resumming the collinear 
logarithms, represented by these additional terms, are modest for this case - a non-obvious 
result on account of the large K-factor and the significant /i-dependence of the FFN scheme 
calculations (see next subsection). The net effect of these correction terms is to increase the 
theoretical cross- sect ion. This is encouraging since the NLO FFN scheme result is known to 
be systematically smaller than the experimentally measured cross section at the Tevatron. 



However, this increase appears to fall short of the current observed discrepancy [0,0]. Cf. 
Fig. |. 

A comparison of Figs. || and shows that, interestingly, the HEl (cthei) contribution 
to the heavy quark production cross section in the heavy quark parton picture is quite 



^^For the purposes of this comparison, we use the same 5-fiavor PDF's for both Fig. ^,b and 
Fig. ^,b so that the only difference is how we combine the terms. Using 4-fiavor PDF's (as 
would be appropriate for the FFN scheme, without the subtraction terms) yields virtually 
indistinguishable curves differing by ~ 1% at pt = 10 GeV, and ~ 3% at pt = 100 GeV. 
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FIG. 8. Comparison with experimental b-production data at the Tevatron = 1800 GeV, 
\y\ < 1. [p6| , p^ The dashed hnes represent the /x- variation in the FFN scheme, and sohd hnes 
represent the //-variation in the VFN scheme. To gauge the // scale variation, we choose ^ = Mt/2 
for upper curves, and // = 2Mt for lower curves. 

comparable to the HCl contribution in the complementary FFN scheme view {cthci, no 
mass-subtraction) - to within about 10%. Thus, at least for this energy range, the heavy 
quark parton picture overlaps considerably with the FFN heavy flavor creation picture as far 
as the inclusive pt distribution is concernedly — these two pictures are complementary rather 
than mutually exclusive, as sometimes perceived in the literature. It is, of course, much easier 
to calculate the tree-level HEl cross section (a text book case) than the HCl one (a tour de 
force). Thus, for this physical quantity, the heavy quark parton picture represents a much 
more efficient way to arrive at the right answer. This approximate equivalence between the 
HEl and HCl contributions to the inclusive pt cross section cannot, of course, be taken 
literally.[^ The two contributions do not have the same (s, pt,fi,y) dependence - in fact, 
the n dependence can be rather different, as we will discuss next. 



^^An equivalent way of saying this is: the subtraction terms, which represent the overlap 
between the two, are a reasonable approximation to both in this energy range. Thus the 
"correction" to either one, represented by the combination of the other with the correspond- 
ing subtraction, are relatively small — as demonstrated above. 

^^The HCl diagram, of course, contains a lot more information on detailed differential 
distributions (such as non-back-to-back-jets) which is not contained in the HEl diagram. 



pp ^bX,Vs= l.STeV, ly^kl 
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B. Scale dependence of the cross section 



In Fig. ^ we show a representative plot of a{pt,^) vs. jj, ai pt = 20 GeV. The tree-level 
HCO and HEl terms give the dominant contributions. In addition to the common 
factor, HCO is predominantly driven by the gluon distribution, and this is a decreasing 
function of yU. On the other hand, the tree-level HEl term is driven by the heavy-quark 
distribution, and this is a increasing function of /i. These two components compensate 
each other. Thus the full LO cross-section in the VFN ACOT scheme has a moderate n 
dependence. 
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FIG. 9. Scale (^) dependence of the leading (a^) order contributions to the cross section 
da /dpi /dy\y=o [nb GeV"^) for b production at 1800 GeV with pt = 20 GeV in the ACOT formalism. 



The situation is different for the FFN scheme, shown in Fig. |T0[ Here one finds that both 
the LO-FFN and NLO-FFN results (proportional to light parton distributions) are decreas- 
ing functions of /i, resulting in a steep /i dependence for the combined result, TOT-FFN. 
If we were to compute higher order corrections in the FFN scheme, we would eventually 
observe compensating terms to reduce the /i dependence (as we know the "all-orders" result 
must be independent of /i). The corrections to the FFN scheme result which have been re- 
summed in the ACOT scheme into HE and GF contributions (minus subtractions) represent 
a part of these higher-order effect. 



In Fig. |TT]we compare directly the TOT-FFN and the TOT- VFN results. The TOT- VFN 
result ranges from ~ 5% to ~ 25% above the TOT-FFN result depending on the choice of /i. 
The additional resummed terms (labeled Difference) is seen to improve the ii dependence. 
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FIG. 10. Scale (/i) dependence of the main contributions to the cross section 

da / dpi / dy\y=Q (n6GeV~^) for b production at 1800 GeV with pt = 20 GeV organized in the 
FFN scheme. 

(See also Fig. ^ below for LHC energies.) Perhaps the improvement is not as complete as 
one might expect. This suggests there is still some non- negligible physics missing from the 
calculation. One possible source might be the NLO-HE process (HE2), on account of the 
large size of HEl. This is an 0{al) contribution which has yet to be computed. 
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FIG. 11. Comparison of the total cross section da / dpj / dy\y=Q (n6GeV~ ) in the ACOT and 
FFN formalism vs. /i for b production at 1800 GeV with pt = 20 GeV. The Difference curve 
represents the additional resummed contributions included in the ACOT result. 



Fig. displays the scaled cross section vs. the physical variable pt for the range of ix 
[Mt/^, 2Mt\. The upper band corresponds to /i = Mt/I, and the lower to /i = 2Mt. We see 
the increase in central value of the cross-section as well as the reduction in scale dependence 
in the ACOT scheme compared to the FFN scheme. The improvement is not dramatic in 



either case at this energy. One expects this to change at higher energies. Fig. |1^ shows the 
corresponding results at the LHC energy of y/s = 14 TeV. For comparison, we also present 



results for b-production at the CERN SppS energy of = 630 GeV in Fig. ^ 
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FIG. 12. Variation of the total cross section da / dpi / dy\y=Q (nftGeV^) in the ACOT and 
FFN formalism vs. pt for h production at 1800 GeV. To gauge the // scale variation, we choose 
^ = Mt/^ for upper curves, and ^ = 2Mt for lower curves. 
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FIG. 13. Same as Fig. 12, with b production at 14 TeV. 
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C. The Rapidity Distribution 



Although the inclusive HEl contribution in the ACOT was roughly comparable to the 
HCl term in the FFN schemes, the dependence on the individual variables (s, pt,fi,y) can 
be quite different. Having already investigated the //-dependence, we turn to the rapidity 
dependence of the underlying processes. 

In Fig. 13 we compare the rapidity distribution for the HEl and the HC0 processes. 
To more easily compare the relative shape, we have scaled the two curves to equal area. 
We observe that the HEl process yields a broader rapidity distribution than the HC0 
processes. In part, this is expected as the HEl process includes a t-channel gluon exchange 
which can give an enhanced contribution in the forward direction. To see the relative effect 
on the rapidity distribution for the complete next-to-leading order calculations, in Fig. ^ 
we display the ratio of the TOT-VFN compared to the TOT-FFN cross section as a function 
of the rapidity. 




FIG. 15. Comparison of the rapidity (y) distribution for the HC0 and HEl processes for b 
production at 1800 GeV with pt = 20 GeV and // = Mt- The curves are scaled to equal area to 
faciUtate comparison of the shapes. 



D. Comment on related work 

It is worth mentioning that the resummation of large logarithms in the fixed-order cal- 
culations into parton distributions and fragmentation functions has also been studied by 
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FIG. 16. Comparison of the total cross section da/dp^/dy in the ACOT and FFN formalism 
vs. y for b production at 1800 GeV with pt = 20 GeV and fi = Mt- The curves are scaled by the 
total FFN cross section to facilitate comparison of the relative magnitude. 



Cacciari and Greco p4[ . In the notation of Eq. 0, their approach amounts to adopting the 
following ansatz for the hadro-production cross section: 



^"ab = E /{?]} ® KUmn = 0) ® rff 

a,b,c 



(17) 



where (a, b, c) are summed over all parton flavors including H, and a'^^^{mH = 0) is given 
by the zero-mass (i.e. light-parton) NLO jet cross section calculation. This is a good ap- 
proximation in the asymptotic region pt ^ but it does not reproduce the right physics 
when Pt is not significantly larger than the quark mass. The main result of their calculation 
was that the predicted 6-production cross section has less scale dependence than the FFN 
scheme result, but it lies within the uncertainty band of the latter. Hence their predictions 
lie substantially below the experimental measurement and somewhat below our results. Our 
theory resums the same large logarithms associated with final state collinear singularities. 
But in addition, (i) we treat initial state parton distributions and final state fragmentations 
of the heavy quark symmetrically; and, more importantly, (ii) by keeping the heavy quark 
mass in the hard cross section according to the general factorization theorem, our results are 
applicable over the entire energy range. Strictly speaking, by using the order jet cross 



section in Eq. Ref. includes higher order corrections to the corresponding terms in 
our calculation (GFl and HFl terms in Eq. However, as shown above, the contribution 
from the latter is already small. Hence, the corrections are not expected to be important. 
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and they are, in any case, of the same numerical order as NNLO flavor creation terms which 
are not currently calculated. Nonetheless, this difference may be responsible for the seem- 
ingly smaller fi dependence of their results. It is hard to be sure because their calculation is 
applicable only when px ^ itlh- 



V. CONCLUDING REMARKS 



In this paper, we have systematically developed the theory for hadro-production of heavy 



quarks according to the natural pQCD scheme of Ref. |jT3|,|T5|,0,|T9[ which generalizes the 
conventional (zero-mass) "improved QCD parton model" to include quark mass effects. 
This formalism has the advantage that it contains the correct physics over a wide range of 
energies and Q (pr), and it reduces to the conventional results both at the low and the high 
energy limits. In particular, it coincides with the widely used FFN scheme in its natural 
region of applicability where Q is or the same order of magnitude as rriH - the one-large-scale 
region. Improvement over the FFN scheme results become important when Q ^ rriH- This 
is manifested both in the reduced theoretical uncertainty and in the increased cross-section. 
The improvement comes at the price of somewhat more complicated calculations. As seen 
in Eq.|l^ and Fig.|, in addition to the FFN scheme contributions, one needs to compute the 
other terms involving flavor excitation and fragmentation processes and their subtractions. 
As emphasized in Sec. |nT| , these calculations must be done with care, due to the delicate 
cancellations required. 

This more complete and consistent theory is, however, not a cure-for-all. Within our 
scheme, the residue scale dependence seen in Sec.|^ may suggest that certain non-negligible 
higher order terms still need to be included. Our cross-section predictions, although higher 
than the existing FFN ones, still fall somewhat short of the b cross-section measured at the 
Tevatron. Some physics effects not included in our formalism could be important: notably, 
those related to large logarithms of the type ln{Q'^/s),ln{mjj/s) ~ Ina; which needs to be 
separately resummed. This is another example of the small-x problem. In this paper, we 
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also have not addressed questions concerning the hadronization of the heavy quark which is 
not yet fully understood. 

As an important physical process involving the interplay of several large scales, heavy 
quark production poses a significant challenge for further development of QCD theory. 
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